AC-field-controlled localization-delocalization transition in one dimensional disordered 

system 



O 
O 

(N 



(N 

CZ5 



Wei Zhang and Sergio E. Ulloa 
Department of Physics and Astronomy, and Nanoscale and Quantum 
Phenomena Institute, Ohio University, Athens, OH 45701-2979 
(Dated: February 6, 2008) 

Based on the random dimer model, we study correlated disorder in a one dimensional system 
driven by a strong AC field. As the correlations in a random system may generate extended states 
and enhance transport in DC fields, we explore the role that AC fields have on these properties. 
We find that similar to ordered structures, AC fields renormalize the effective hopping constant to 
a smaller value, and thus help to localize a state. We find that AC fields control then a localization- 
delocalization transition in a given one dimensional systems with correlated disorder. The compe- 
tition between band renormalization (band collapse/dynamic localization), Anderson localization, 
and the structure correlation is shown to result in interesting transport properties. 
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I. INTRODUCTION 

The dynamics of electrons in ordered semiconduc- 
tor superlattices driven by electric fields has received a 
great deal of attention. For example, well-known predic- 
tions for quantum mechanical behavior, such as Wannier- 
stark ladders and Bloch oscillations, which are difficult 
to observe in ordinary solids, were verified in beauti- 
ful experiments. 1 Furthermore, behavior such as nega- 
tive differential conductance^ fractional Wannier-Stark 
ladders,- excitonic Franz-Keldysh effect^ among others, 
have been the focus of recent work. Of particular re- 
cent interest are the localization and derealization be- 
havior of electrons in the presence of external electric 
fields, which have direct effect on the macroscopic trans- 
port properties of the system. For example, in a system 
with both AC and DC fields, an appropriate AC field 
will delocalize the Wannier-Stark ladder states induced 
by a strong pure DC field. 5 An intense AC field can by 
itself also lead to dynamical localization of the carriers. 6 
This phenomenon has been studied extensively in vari- 
ous systems^ especially in quantum dot pairsf 8 and finite 
linear arrays^ 

Certainly, disorder or imperfections are unavoidable in 
a real system. It is widely known that localization due to 
disorder plays a fundamental role in a variety of physical 
situations. In particular, it has been of interest to investi- 
gate disordered systems in the presence of electric fields: 
Hone et alm^ and Zhang et a/ii studied the case of one 
impurity in the presence of AC fields. Holthaus et al. 
studied AC-field-controlled Anderson localization in dis- 
ordered semiconductor superlatticesii^ As scaling theory 
shows that all eigenstates in disordered one-dimensional 
(ID) systems are localized^ previous studies have fo- 
cused on the effects of electric fields on the localization 
length iiS 

On a related area, the existence of metallic states in 
a class of conducting polymers, such as polyaniline and 
heavily doped polyacetylene, was identified by Dunlap 



et alM> with extended states in ID systems that exist if 
short-range correlations in the disordered structure are 
taken into account. The existence of extended states in 
this random dimer model was also verified in experiments 
with GaAs-AlGaAs superlattices designed to exhibit such 
correlated disorderii 5 . 

More general correlations have also been studied in 
ID systems. Perturbation theories for the random dimer 
model were developed in Ref. phi ], and particle trans- 
port in models with correlated diagonal and off-diagonal 
disorder were discussed also by Floresiii The random 
dimer model driven by a DC field was studied ini&. The 
derealization behavior in ID models with long-ranged 
correlated disorder for the on-site energies was stud- 
ied by a renormalization technique^ and by a Hamil- 
tonian approach^ More recently, the Kronig-Penney 
model with correlated disorder was studied^ demon- 
strating that a mobility edge may exist for disordered 
systems with appropriate long-range correlated disorder. 
The role of structural correlations in the sequence disor- 
der in DNA molecules has also been studied recentlyj22i2ii 
as this represents a real system that exhibits structural 
correlations between the diagonal and off-diagonal ele- 
ments in a tight-binding representation. 

In this rich context, it is important to study the com- 
petition between dynamic localization (due to AC field), 
Anderson localization (due to disorder) and the correla- 
tion in the disorder, and to investigate the role of ex- 
ternal fields on the transition from the localized to the 
delocalized state in ID systems. In this paper, we con- 
centrate on the short range correlations of the random 
dimer model, and study the localization-delocalization 
transition driven by external AC electric fields. We find 
that AC electric fields induce a transition from extended 
to localized states under suitable conditions, and find the 
transition point analytically in the high frequency limit. 
We also show that the transition for lower frequencies is 
shifted in field, as a precursor of DC-field results. Al- 
though our results are for a relatively simple ID model 
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potential, we expect that they will be relevant for a va- 
riety of dissimilar systems, including polymers^ exciton 
transfer in active media^i semiconductor superlattices^ 
quantum dot arrays and even hole transport in complex 
molecules 2 ^^, as their dynamics is described well by ef- 
fective ID models. 

After introducing our general model in Sec. II and pre- 
senting an analysis of the high frequency regime for the 

I 

A. Model 

We consider a ID random dimer model driven by an AC electric field. The appropriate Hamiltonian is then 

H = Em{ £ ™ a m a m + R ( a m+i a m + a m a+ n+1 ) + me<iE(*)a+ a m }, (1) 



smgle-dimer system m Sec. Ill, we present numerical re- 
sults and general discussions in Sec. IV. 



II. RANDOM DIMER MODEL IN AC 
ELECTRIC FIELDS 



where R is the hopping amplitude between nearest neigh- 
bors, E(t) — E\ cos(wi) is the time-dependent field with 
frequency u>, d is the constant separation between chain 
sites, and the on-site energy parameter is e m = e a (or £{,) 
with probability Q (or 1 — Q; here we typically choose 
Q = 1/2, as it represents the most disordered system, 
although other values are also used), and Sb is assigned 
to a pair of nearest neighbor sites when it occurs. Since 
the Hamiltonian is periodic in time, the Floquet theorem 
implies that the state can be written as 

^M) = e- fe *£„C n (i)<^(x), (2) 



where e is the quasi-energy, 4> n is the Wannier state and 
C n is the probability amplitude for an electron on site n 
at time t, which is periodic in time, i.e., C n (t) = C n (t + 
T), with T = 2tt/uj. The Schrodinger equation 



= W{x,t) (3) 



can then be written as 



ij%C n (t) = (e„ - e)C n {t) + R(C n+1 + C n _ t ) + nedE 1 cos(ojt)C n . (4) 

Since the term containing the electric field is proportional to n, it is not suitable to perform perturbative calculations. 
We introduce the following transformation 

C' n (t) = C n {t)e inl3s ' ln ^ , (5) 
where (3 — edEi/uo. It is easy to see that |C4(t)| = |C„(t)|, C' n (t + T) = C' n (t), and that C' n (t) satisfies the equation 

if t C' n {t) = (e n - e)C'Jt) + R( e -^^C' n+1 + e^WC^). (6) 
Since C' n (t) is periodic in time, we can expand it in Fourier series, C' n (t) — J2 m A™e lmult . Using the identity 25 

e « s in W = ^+~_ oo (_l)"» J m (z)e"™ 9 , (7) 
where J rn is the m-th Bessel function, we can obtain an equation for A™ 

(Sn-e + m^ + R^ [AZ+K-lYM-fi+A^K-iyMP)] =0. (8) 

B. High frequency limit 

We analyze here the behavior of our model in the high frequency regime. Apart from illustrative, this limit is of 
practical importance, since much interest exists in systems driven by THz fields. 4,26 In the high frequency limit, 
gives the most important contribution to C n (t). The equation for is 

(e n - e)Al +RJ2 l (A l n+1 (-l) l J l (-P) + A^-l) 1 J0)) = 0. (9) 
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We first keep only the terms with 1 = 0, and obtain, 

{e n -e)Al = RJ (/3)(A° +1 + Al_ 1 ). 



(10) 



This equation indicates simply that in the high frequency limit, the effect of the AC field is to suppress hopping and 
change R to an effective hopping constant R e g = RJq{(3). This is a well-known result in driven systems^ 
Let us consider the corrections coming from terms containing A^ 1 - The relevant equations are 



(e„ - e)Al + RJ (p)(A° n+1 + A t) n _ 1 ) + RJ^Al 



a; 



A 1 



a: 



o 



(e n - e + lu)AI + RJ (p)(Al +1 + A^) - RJ^A* - A^) = 



and 



(e n - e - u)A-i + RJomA-l, + A~\) + R.h{(3)(Al +1 - A^) = . 
These equations (|11H can be further simplified, and in conjunction with (|10|) . we find that 



(B n+ i--B«-i)(l-(^ 



+ 



(RJo 



n+3 



B n - 3 ) 



2R JpJ\ 



2R JqJi 



(^n+3 + A-n-z) 



(11) 



(12) 



r 



where we have defined B n = A x n — A~ . Thus the higher 
Fourier component corrections to the effective bandwidth 
RcS = RJq([3) are of higher order in R/ui, which of course 
makes them small in the high frequency limit, R/oj ^ 1. 
In this limit, the model with AC field behaves essentially 
as a system without electric field, except for a rescaling 
of the bandwidth given to the lowest order by i? c ff = 

In the random dimer model without AC field, the 
localization-delocalization transition occurs when e_ = 
\e a — Sb\ is twice the bandwidthiii We then intuitively 
expect that e_ = 2i? c ff would be the transition point be- 
tween localized and delocalized states in the presence of 
the AC field. Before giving numerical evidence for this 
transition in the full random dimer system, we analyze 
the single impurity case to gain further understanding of 
this problem. 



III. SINGLE IMPURITY-DIMER CASE 

It is instructive to see what happens when only one 
impurity-dimcr is involved in an otherwise periodic ID 
chain. We consider both the cases with and without an 
AC field for comparison. 



A. Static case 

Let us first consider the scattering effects introduced 
by a single site-impurity in a chain, in the absence of AC 
fields We let all site energies be e a , except for site 



where it is Sb- From the eigenvalue equation 

(e n — e)C n + R(C n+ i + C„_i) = 0, 
we can get the transmission probability 



IT 



2 _ (4R sin kd) 2 
~ e 2 _ + (ARsinkd) 



7 ' 



(13) 



(14) 



where k is the wave vector of the incoming Bloch wave, 
and e_ = |efc— £ a | measures the impurity detuning, and as 
such is a measure of the "disorder" or impurity strength. 
One can see that for this one-impurity case, \T\ 2 < 1 for 
£- =/= 0. In a random multi-impurity system, a series of n 
scattering events would naturally lead to \T\ 2n 1, and 
very small amplitude for the outgoing wave, resulting in 
localization in the thermodynamic limit. 

If instead, we assign a pair of sites and 1 with en- 
ergy 6b to form a single dimer impurity, the transmission 
probability isi^ 

e 2 _{e- + 2R cos kd) 2 



TV = l- 



e_{e- + 2Rcos kd) x 



4i? 2 sin^ kd 



(15) 



In this one impurity- dimer case, |T| = 1, whenever 
e_ = —2R cos kd. This is a sort of resonance effect due 
to the internal structure of the impurities. Thus, in the 
presence of the peculiar kind of short-range correlated 
disorder described by the random dimer model, there are 
states with unity transmission probability, which clearly 
have an extended character (even if they only appear at 
a single value of the energy). This one impurity-dimcr 
calculation, makes intuitive the appearance of extended 
states in the random multi-dimer system at peculiar en- 
ergy values, and captures qualitatively the reason for the 
unusual behavior of the random dimer model Il4ll27l. 
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B. The case with AC electric field 

We now turn to explore what happens when an AC 
electric field is turned on. Hone et al. studied the system 
of an isolated defect driven by strong electric fields^ We 
will make use of Green's functions in terms the Floquet 
formalism. We consider the resolvent operator as a func- 
tion of the complex frequency z 



G(z) 



1 

z - H 



G°{z) + G°{z)VG{z), 



(16) 



where H = Ho + V, Hq is the unperturbed Hamiltonian, 
V is the impurity potential, and G°(z) 



1 



rr- ■ In the 

Z - Hq 

representation of Wannier states, the Green's function is 



G jl = G% + Y,G%,V rv G n . 
In the high frequency limit, R/ui <C 1, G° is^i 



nu - _J_£ 

^ ~ i?cff 1-q 2 ' 

where R c ff = RJq ( 6C [fi j , and 



(17) 



(18) 



2i? cf f 



(19) 



The sign is chosen such that q falls inside the unit circle. 
For an isolated defect with Vji = v5j^8o,i (where v = 
E- = \eb — £ a \ is the case of different site energy), the 
probability p(l) for the defect state (with energy e — 
R c s(q + i)) to occupy the ^ th site is determined by the 
residue of Gu . One finds thafei^ 



p(0 



4R 2 



AR, 



1 - 



eir 



2R 



elf 



2\l\ 



(20) 

Notice that p(l) falls exponentially from the site I = 
where the defect is localized, with a characteristic de- 
cay length that is reduced for increasing \v\/R e s, as one 
might suspect. 

Now we consider a one-dimer model, with Vji = 
v(6joSoi + SjiSu), and v = e_. After some calculation, 
we find 



Gu = G° u + vG a w Goi + uG^Gu 



Gm = 



(1 - vG^Gl 
tD 



7^0 
^01 



(l-uGp(l-vG^)- Q ^ 
(-• _ vGiqGqx + (1 - vGqq)G 11 



(21) 



so that 



G u = aq+ ^va 2 q 2L 



(1 + g) 2 



(1-g) 2 



1 - aq(l + q) 1 — aq(l — q) 



(22) 



where a is defined as 



for the pole at qi — 



1 



27TT' we § et 



(23) 



and since q\ < 1, this corresponds to a localized state. 



For the pole q2 = 



2t~t > one e ets 
p(O-fe) 2 '. 



(24) 



When ^ = e_ > 2R c g, q^ < 1, and this corresponds to 
a localized state with localization length ~ 1/ ln(27 — 1). 
When 7 approaches 1 from above, the localization length 
diverges, indicating a transition to a delocalized state. 
Thus the single dimer impurity in an AC field yields the 
same conclusion as in the high frequency case, and the 
transition from localized to extended states occurs at the 
point v — 2i? c ff. 



IV. NUMERICAL RESULTS AND DISCUSSION 

Most of our numerical calculations were performed on 
a chain with 1501 sites. We solved equation 10} with 
initial condition C n (t = 0) = <5„.o, and analyze the sub- 
sequent development. The site energies were chosen from 
a bi- valued distribution, e n = e a and e n — Sf,, with prob- 
ability |. As the site prob remains, the "degree of dis- 
order" is controlled by larger values of £_ = |e& — e a \, as 
we will see in what follows. 



A. High frequency regime ^ < 1 

In Fig. 1 we show numerical calculations of the mean- 
square displacement, (m 2 ) = J2 n n \C n \ , versus time. 
One can see in curve (a) that when e_ = |e& — e a \ = R e g, 
the mean-square displacement (m 2 ) ~ i 3 / 2 . This is 
known as the super diffusive transport regime. Diffu- 
sive transport {{m 2 ) ~ t) is shown either when e_ = 
0.97 • 2i? c ff (curve b) or perhaps for e_ = 2i? e ff (curve 
c). On the other hand, curve (d) for e_ = 2R = 2.27i? c ff 
shows how the mean square displacement is increasingly 
bounded (subdiffusive), (m 2 ) ~ t 36 . Further increasing 
e_ > 2R e s, as shown in curve (e)(e_ = u = 3.78i? G ff) 
results in completely bounded motion (m 2 ) «const., as 
anticipated from the analytical discussion above. We be- 
lieve the subdiffusive behavior for e_ = 2R > 2R e g is a 
crossover behavior due to finite size effects, masking the 
anticipated extended — > localization transition at e_ = 
2R e g. In fact, from the Eq. (15) (with R replaced by 
R e g) for the transmission probability, we may estimate 

the localization length A ~ -1/ In |T(/c)| 2 ~ - — ^r?, 

where e_ = —5 — 2R c g, and S is small. The number 
of "extended states" (states with localization length A 
larger than the system size L) is AJV ~ L(-j- — 5) 1 / 2 . 
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(a) 



Dotted reference line: slope=-1/2 




FIG. 1: The mean-square displacement divided by (Rt) 3 ^ 2 
for varying amounts of disorder in the high frequency regime 
R/lo — 0.3. Amplitude of AC electric field is j3 = edEi/cu — 
0.7; effective hopping constant R c g = 7?Jo(/3) = 0.88127?. (a) 
With e_ = \sb — e a \ = Reft results in superdiffusive transport, 
(m 2 ) ~ t 3/2 ; while (b) = 0.97 • 27? cff , and (c) e_ = 27? cff 
result in diffusive behavior, (m 2 ) ~ t, compared with dotted 
line of slope -1/2. For (d) e_ = 27?, we see (m 2 ) is nearly 
bounded, and for (e) e_ = u = 3.787? e gf, we find completely 
bounded {m 2 }, compared with dashed slop -3/2. 



We believe theses "extended states" lead to the subdif- 
fusive behavior. With increasing of S, as in curve (e), 
or increasing system size, one gets AN = 0, if -h < S. 
It is very difficult to go beyond the crossover regime by 
numerical simulations, as it requires simulations in very 
large system sizes, with longer equilibration times, re- 
quiring longer simulation time to obtain accurate values 
of the self-averaged quantities in a In Rt fashion. Scaling 
studies of this transition would be interesting. 

The structure of Fig. 1 is similar to that shown in 0] 
in the absence of AC fields. For a fixed AC electric field 
amplitude, there is a transition from extended to local- 
ized state behavior with increasing disorder. The role 
of AC electric fields can be seen to effectively decrease 
the hopping constant, thus contributing to the localiza- 
tion of carriers. For example, for e_ < 27?, in the case 
without electric fields, results in extended states (and dif- 
fusive transport)^ However, when an AC electric field 
(with edEi/uj — 0.7) is turned on, the mean squared dis- 
placement is suppressed. This localization-delocalization 
transition is clearly induced by the AC electric field, as 
the transition shifts to e ~ 2i? e g 

It is interesting to see the situation for a stronger field 
(3 = edEi/uj = 2.405, (the first root of J ). In this case, 
Reft — 0, and we expect that even for very weak dis- 
order the state will be localized. Our results in Fig. 2 
(with very small disorder, e_/7? = 0.33) show that this 
is indeed the case. This is also in agreement with the fact 



that even in the limit of e_ = 0, i.e. when there is no dis- 
order, the states are localized when band collapse occurs 
(i.e. 7? c ff = 0). 28 This is nothing but the well-known dy- 
namical localization. 6 In Fig. 2 one notices that there are 
oscillations in the mean square displacement. This is the 
manifestation of the time dependence of the electric field 
in this case of weak disorder. In fact these oscillations 
also exist in Figs. 1 and 3, except that they are nearly 
invisible in those cases because the disorder e_ and the 
displacements being larger, "hide" the oscillations. 




FIG. 2: The mean-square displacement divided by (Rt) 3 ^ 2 for 
P — edEi/uo = 2.405, Jo(/3) = 0, in the high frequency regime, 
R/lo = 0.3. The effective hopping constant is 7? c ff = 0, and 
e~ = \sb — £a\ = R/3. Notice even weak disorder results in 
bounded displacement, characteristic of localization. 



B. Low frequency regime ~ 1 



In Fig. 3 we show the transition from localized to ex- 
tended state in the low frequency limit. One can see that 
the transition point is no longer the same as in the high 
frequency regime e_ = 2RJo(edE\ / 'lj) . For example, for 
£_ = 7? ff < 2R c g (curve b) the state is sub-diffusive, 
with (m 2 ) ~ 7 67 . As expected, when e_ is small enough, 
for example e_ = 0.347? c gf (curve a), the state is extended 
and shows superdiffusive behavior. These results are of 
course different from the high frequency limits, since in 
this regime our previous analysis fails. Furthermore, in 
the extreme low frequency limit, the system tends to that 
with a DC field: It is known that localized states have a 
power-law behavior in DC field, instead of the more typ- 
ical exponential localization in ID disordered systems^ 
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In(Rt) 



FIG. 3: The mean-square displacement divided by (Rt) 3 ^ 2 
for varying amounts of disorder in the low frequency regime 
R/lu = 1. The amplitude of AC electric field is edEi/uj — 
0.7. The effective hopping constant is R c g = 0.8812i?. (a) 
£_ = 0.34R c ff shows extended, superdiffusive behavior; (b) 
£_ = Rtfi shows near localization well below the critical value 
of e_ = 2R c g for high frequency. 



C. Inverse participation ratio 

To understand more clearly how the AC electric field 
controls the degree of localization, it is useful to extract 
information from the Floquet states for the system driven 
by periodic electric fields. The Floquet states be 
expanded with respect to the Wannier states <pi, 



(25) 



We calculate the averaged inverse participation ratio P, 



N 



1 Jo dt \ 



,{m) 



(*)l 



(26) 



If a Floquet state is nearly localized at individual Wan- 
nier states, P tends to 1, while P vanishes as 1/N if 
the state is extended; the larger P characterizes a more 
localized state. In Fig. 4, we show P for different val- 
ues of £_ = \eb — s a \ and dimer concentrations versus 
electric field strength E\ in the high frequency regime, 
R/lu = 0.1. We find sharp peaks at edEi/uj — 2.405, 
as this value results in R e g — 0, and thus the effective 
hopping along the chain vanishes. We can enhance the 
degree of localization in the random dimer model by in- 
creasing the detuning e_ or the dimer concentration Q. 
For cases (a) and (b) in Fig. 4, Q is the same (= 0.5), but 
e_ changes from 0.16 in (a) to 0.07 in (b). In contrast, 
for (b) and (c), the value e_ is the same, but Q = 0.2 is 
smaller in (c). It is clear that P is larger overall for the 
more disordered systems, and although a peak appears 
always at edEi/uj ~ 2.4, decreasing disorder suppresses 
the peak value and overall amplitude of P. One can also 



observe that there is a relatively sudden enhancement of 
P for the system in (a) for edEi/u > 0.9, while for (b) 
and (c) this occurs between edEi/ui ~ 1.7, and 3.3. From 
our previous discussion, we know that the localization- 
delocalization transition occurs at e_ = 2RJo(edEi/uj). 
From this formula, we find that the transition point for 
(a) is in fact at edE\juj = 0.92, while for (b) and (c) it 
occurs for edEi/u = 1.78, and 3.33. These match very 
well with our numerical calculation. 
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FIG. 4: The averaged inverse participation ratio versus elec- 
tric fields for the random dimer model. Lattice size N — 81 
and R/lu — 0.1. (a) e_ = 0.16, the concentration of dimer 
of energy e b is 50% (Q = 0.5); (b) e_ = 0.07 and Q = 0.5; 
(c) e_ = 0.07, but with concentration of dimer Q = 0.2. No- 
tice localized peak at edEi/uj = 2.4 has decreasing amplitude 
with decreasing e_ or Q. Region of high P values agrees with 
expected e_ = 2RJo(edEi/u>) (see text). 



To elucidate further the role of correlations, we com- 
pare P in an Anderson model (without correlations) with 
a random dimer model system, as shown in Fig. 5. For 
a more quantitative comparison, we let the variance of 
the Anderson model distribution, W 2 /12, be the same 

as that in the random dimer model, ^(e^ + ef). It is 
evident that P is much larger in the Anderson model 
(indicating a more localized system), and that P varies 
smoothly with electric field, indicating no localization- 
delocalization transition with fielded This figure also in- 
dicates that an important effect of the presence of the 
random dimer short range correlations is to delocalize a 
few states, reducing globally the value of P in the sys- 
tem. It is clear that the dynamical behavior of a system 
with correlated disorder is a subtle competition between 
correlation and disorder. 
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FIG. 5: The averaged inverse participation ratio versus elec- 
tric fields for random dimer model. Lattice size N — 81 
and R/uj — 0.1. (a) For Anderson model with W = 
y6(£a+ej;) = 0.392; (b) for random dimer model with 
£b — £a = 0.16, the concentration of dimer of energy et is 
50%. It is clear that the Anderson model with uncorrelated 
disorder is more localized and has larger P values at all fields. 



D. conclusions 

We have studied the AC-field controlled random dimer 
model. The dynamics of our system depends on the 



competition between band renormalization (band col- 
lapse/dynamic localization), Anderson localization, and 
the correlation (dimer structure). We find that there is 
an AC electric field induced transition from extended to 
localized states, which is absent in the Anderson model. 
The transition point is found analytically for the high fre- 
quency limit, and found to occur when e_ = \sb — s a \ ~ 
i? ff = 2RJa(edEi/ lo). The dynamical localization is not 
only recovered as a natural limit in the absence of disor- 
der, but also shows its effects in the transport properties 
of the system with disorder and correlation (the peaks 
in Figure 4 and 5). The generalization of our results 
to a V-dimer model is straightforward, and expected to 
yield qualitatively similar results. Our theoretical pre- 
dictions could be checked in a variety of systems, and es- 
pecially on experiments in GaAs-AlGaAs random-dimer 
superlattices^ In experiments, tuning external AC field 
is a relative easy task compared with changing disorder 
or correlation in a desired way. Generalizations to differ- 
ent and more complex correlations are also expected to 
give interesting results. 



Acknowledgments 

We acknowledge helpful discussions with J. M. Villas- 
Boas, and support from the US-DOE and the 21 s * Cen- 
tury Indiana Fund. 



1 E. E. Mendez, F. Agullo-Rueda, and J. M. Hong, Phys. 
Rev. Lett. 60, 2426 (1988); E. E. Mendez and G. Bastard, 
Phys. Today 46, 34 (1993); C. Martijin de Sterke, J. N. 
Bright, P. A. Krug, and T. E. Hammon, Phys. Rev. E 
57, 2365 (1998); S. R. Wilkinson, C. F. Bharucha, K. W. 
Madison, Q. Niu, and M. G. Raizen, Phys. Rev. Lett. 76, 
4512 (1996); M. B. Dahan, E. Peik, J. Reichel, Y. Castin, 
and C. Salomon, Phys. Rev. Lett. 76, 4508 (1996). 

2 R. Tsu and G. H. Dohler, Phys. Rev. B 12, 680 (1975). 

3 X.-G. Zhao, R. Jahnke, and Q. Niu, Phys. Lett. A 202, 
297 (1995). 

4 K. B. Nordstrom, K. Johnsen, S. J. Allen, A.-P. Jauho, 
B. Birnir, J. Kono, T. Noda, H. Akiyama, and H.Sakaki, 
Phys. Rev. Lett. 81, 457 (1998). 

5 M. Holthaus and D. W. Hone, Phil. Mag. B 74, 105 (1996). 

6 D. H. Dunlap and V. M. Kenkre, Phys. Rev. B. 34, 3625 
(1986); Phys. Lett. A 127, 438 (1998); X.-G. Zhao, Phys. 
Lett. A 155, 299 (1991); 167, 291 (1992). 

7 See, for example, M. Grifoni and P. Hanggi, Phys. Rep. 
304, 229 (1998). 

8 J.M. Villas-Boas, W. Zhang, S.E. Ulloa, P.H. Rivera, and 
N. Studart, Phys. Rev. B 66, 085325 (2002). 

9 J.M. Villas-Boas, S.E. Ulloa, and N. Studart, Phys. Rev. 
B 70, 041302 (2004). 

D. W. Hone and M. Holthaus, Phys. Rev. B 48, 15123 
(1993). 



11 A.-Z. Zhang, P. Zhang, S.-Q. Duan, X.-G. Zhao, and J.-Q. 
Liang, Phys. Rev. B. 63, 045319 (2001). 

12 M. Holthaus, G. H. Ristow, and D. W. Hone, Phys. Rev. 
Lett. 75, 3914 (1995). 

13 E. Abrahams, P. W. Anderson, D. C. Licciardello, and T. 
V. Ramakrishnan, Phys. Rev. Lett. 42, 673 (1979). 

14 D. H. Dunlap, H-L. Wu, and P. W. Phillips, Phys. Rev. 
Lett. 65, 88 (1990); P. Phillips and H.-L. Wu, Science 253, 
1805 (1991). 

15 V. Bellani, E. Diez, R. Hey, L. Toni, L. Tarricone, G. B. 
Parravicini, F. Dominguez-Adame, and R. Gomez-Alcala 
Phys. Rev. Lett. 82, 2159 (1999); V. Bellani, E. Diez, A. 
Parisini, L. Tarricone, R. Hey, G. B. Parravicini, and F. 
Dominguez-Adame, Physica E 7, 823 (2000). 

16 J. C. Flores and M. Hilke, J. Phys. A: Math. Gen. 26, 
L1255 (1993); A. Bovier, J. Phys. A: Math. Gen. 25, 1021 
(1992). 

17 J. C. Flores, J. Phys. : Condens. Matter 1, 8471 (1989). 

18 F. Dominguez-Adame, A. Sanchez, and E. Diez, J. Appl. 
Phys. 81, 777 (1997). 

19 F. Moura and M. Lyra, Phys. Rev. Lett. 81, 3735 (1998); 
ibid 84, 199 (2000). 

20 F. M. Izrailev and A. A. Krokhin, Phys. Rev. Lett. 82, 
4062 (1999). 

21 F. M. Izrailev, A. A. Krokhin, and S. E. Ulloa, Phys. Rev. 
B 63, 041102(R) (2001). 



8 



W. Zhang and S. E. Ulloa, Phys. Rev. B 69, 153203 (2004), 
and in preparation. 

D. K. Klotsa, R. A. Romer, and M. S. Turner, Biophys. J. 
89, 2187 (2005), and cond-mat/0508720 

exciton hopping by Dominguez-Adame in PRL ... 
Handbook of Mathematical Functions, M. Abramowitz and 
I. A. Stegun, eds. (Dover, New York, 1965). 

E. H. Cannon, K. N. Alekseev, F. V. Kusmartsev, and D. 
K. Campbell, Europhysics Lett. 56, 842 (2001). 



S. N. Evangelou and A. Z. Wang, Phys. Rev. B 47, 13126 
(1993). 

M. Holthaus, Phys. Rev. Lett. 69, 351 (1992); W. Zhang 
and X.-G. Zhao, Physica E 9, 667 (2001). 
C. M. Soukoulis, J. V. Jose, E. N. Economou, and P. Sheng, 
Phys. Rev. Lett. 50, 754 (1983); F. Delyon, B. Simon, and 
B. Souillard, Phys. Rev. Lett. 52, 2187 (1984). 



